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■ Abstract 

Critical slowing down of the relaxation of the order parameter has phenomeno- 
logical consequences in early universe cosmology and in ultrarelativistic heavy 
ion collisions. We study the relaxation rate of long-wavelength fluctuations of 
XjT^ • the order parameter in an 0{N) scalar theory near the critical point to model 



0^ 



the non-equilibrium dynamics of critical fluctuations near the chiral phase 



I transition. A lowest order perturbative calculation (two loops in the coupling 

I A) reveals the breakdown of perturbation theory for long-wavelength fluctu- 

' ations in the critical region and the emergence of a hierarchy of scales with 

^ • hard q >T, semisoft T ^ q ^ XT and soft AT ^ q loop momenta which are 

^ . widely separated in the weak coupling limit. A non-perturbative resummation 

D i is implemented to leading order in the large limit which reveals the infrared 

' renormalization of the static scattering amplitude and the crossover to an ef- 

fective three dimensional theory for the soft loop momenta near the critical 
^ ■ point. The effective three dimensional coupling is driven to the Wilson-Fisher 

^ . three dimensional fixed point in the soft limit. This resummation provides 

an infrared screening and for critical fluctuations of the order parameter with 
wave- vectors AT ^ k ^ k or near the critical temperature AT ^ mf ^ kus 
with the ultrasoft scale kus = ^ exp [— ^] the relaxation rate is dominated 
by classical semisoft loop momentum leading to r(A;,T) = XT/{2ttN). For 
wavectors k <C k^s the damping rate is dominated by hard loop momenta 
and given by r(A;, T) = 47rT/ [3A^ln(T/A:)]. Analogously, for homogeneous 
fluctuations in the ultracritical region <C k^s the damping rate is given 
by ro(mr,T) = AttT/ [3A^ ln(T/mT)]. Thus critical slowing down emerges for 
ultrasoft fluctuations. In such regime the rate is independent of the coupling 
A and both perturbation theory and the classical approximation within the 
large N limit break down. The strong coupling regime and the shortcomings 
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of the quasiparticle interpretation are discussed. 



lypeset using REVI^ 



I. INTRODUCTION AND MOTIVATION 



The program of relativistic heavy ion colhsions both at Brookhaven and at Cern seeks 
to understand the phase diagram of QCD in conditions of temperatures that were achieved 
during the first 10/is after the Big Bang or densities several times that of nuclear matter which 
could exist at the center of neutron stars. Current theoretical understanding |I| leads to the 
conclusion that QCD could undergo two phase transitions: a confinement-deconfinement (or 
hadronization) and the chiral phase transition. Current lattice data seem to suggest that 
both occur at about the same temperature Tc ~ IGOMeV The consensus emerging in the 
field is that several types of observables will have to be studied simultaneously and event- 
by-event analysis of data will have to be carried out to extract unambiguous signals both 
hadronic and electromagnetic to reveal the presence of a Quark-Gluon Plasma phase. Recent 
results reported from CERN-SPS seem to indicate a strong evidence for the existence of 
the QGP in Pb-Pb collisions, and RHIC at Brookhaven will begin operation soon with Au- 
Au collisions with four dedicated detectors capable of event-by-event analysis of hadronic 
and electromagnetic observables. 

For QCD with only two flavors of massless quarks (u,d) it has been argued that 
the chiral phase transition at finite temperature but vanishing baryon number density is of 
second order and described by the universality class of 0(4) Heisenberg ferromagnets. It 
has also been suggested recently that at finite baryon density there is a second order critical 
point described by the Ising universality class ||5|. Second order critical points are char- 
acterized by strong critical long-wavelength fluctuations and a diverging correlation length 
that could lead to important experimental signatures 0. These signatures would be akin 
to critical opalescence near the critical point in binary fluids and could be observed in 
an event-by-event analysis of the fluctuations of the charged particle transverse momentum 
distribution (mainly pions) 0. These fluctuations are characterized by the typical correla- 
tion length of the order parameter and it has been suggested that the phenomenon of critical 
slowing down, ubiquitous near the critical point of second order phase transitions, can lead 
to strong departures from equilibrium that will determine the value of the correlation length 
when long- wavelength fluctuations freeze out [0] . Critical slowing down of long- wavelength 
fluctuations near a second order critical point is the statement that the long-wavelength 
Fourier components of the order parameter relax very slowly towards equilibrium In 
mean-field theory in classical critical phenomena, the relaxation time diverges proportional 
to the susceptibility near the critical temperature but thermal fluctuations renormalize the 
relaxation time to be of the form T{k = 0) oc .^^ with ^ the correlation length, or at critical 
point T{k) oc with z a dynamical critical exponent [^0. Another similar manifestation of 
an anomalously slow relaxation of long-wavelength fluctuations arises in weakly first order 
phase transitions when the system enters into the mixed phase where the (isothermal) speed 
of sound, which determines the velocity of propagation of long-wavelength pressure waves, 
becomes anomalously small resulting in the softest point of the equation of state. In this case 
there have also been suggestions that there are experimental consequences of this softening 
in relativistic heavy ion collisions in observables related to collective flow and the transverse 
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momentum distributions of particles at freezeout |T^ . 

In classical normal fluids near the critical point the vanishing of the (isothermal) speed of 
sound, critical opalescence (strong scattering of light by long-wavelength fluctuations) and 



critical slowing down are all related [|TT|, and in ferromagnets the spin diffusion constant 
vanishes near the critical point again signaling critical slowing down |T1[ . 

The softening of the equation of state near the critical point of QCD could also have im- 
portant cosmological implications. When the the QGP enters the mixed phase with hadrons, 
the speed of sound becomes anomalously small and the time scale for propagation of pressure 
waves over a given critical wavelength becomes longer than the free fall time for gravitational 
collapse which is then unhindered by the pressure of the hadronic gas. This could lead to the 
formation of primordial black holes [|l^] with a possible imprint in the acoustic peaks in the 
cosmic microwave background |TB|. Other possible cosmological relics from the QCD phase 
transition with a mixed phase had been predicted, from strange quark nuggets to MACHO's 



A familiar argument is typically invoked to state that while the QCD phase transition 
in the Early Universe occurred in local thermodynamic equilibrium (LTE) this may not be 
the case in Relativistic Heavy Ion Collisions. The argument compares the typical collisional 
relaxation time scale obtained from a strong interaction process Tcou ~ 10^^^ sees to the 
time scale for cooling near the critical temperature ~ 160MeV, i.e. T/T ~ H^^ ~ 10~^ sees. 
The argument is that since Tcou ^ H"^ the phase transition occurs in LTE in cosmology, 
whereas in relativistic heavy ion collisions at RHIC and LHC energies these time scales will be 
comparable. However this argument completely neglects the possibility that long-wavelength 
fluctuations could undergo critical slowing down and freeze out, i.e. fall out of local thermal 
equilibrium, even before the phase transition. The freeze-out of long-wavelength fluctuations 
during the phase transition could result in important non-equilibrium effects on the size 
and distribution of primordial black holes or any other cosmological relic just as they could 
lead to important observables in the momentum distributions of charged pions in relativistic 
heavy ion collisions [|6|,|7|,p!0| . 

Indeed there are simpler experimental situations where this is the case, in typical normal 
fluids a collisional relaxation time (away from the critical point) is of the order of 10~^ sees 
while near the critical point (even at 10% of the critical temperature) critical slowing down 
becomes very dramatic and thermalization time scales become of the order of minutes if not 



hours ll,lq,p!i 



Thus the phenomenological importance of critical slowing down for the QCD phase tran- 
sition both in Relativistic Heavy Ion Collisions as well as in Early Universe Cosmology 
motivates us to study this phenomenon in a model Quantum Field Theory that bears on the 
low energy (chiral) phenomenology of QCD, the 0{N) linear sigma model. Furthermore the 
study of critical slowing down is the precursor to a more complete program to understand 
transport phenomena and the relaxation of hydrodynamic modes at or near a critical point. 

Goal: Our goal is to provide a consistent microscopic description of critical slowing 
down at or near criticality directly from an underlying quantum field theory that is at least 
phenomenologically motivated to study the QCD phase transitions. This will be a first step 
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in a program that seeks to offer a consistent description of transport near critical points 
that eventually may be merged with a hydrodynamic description to obtain a more reliable 
picture of critical phenomena near the deconfinement and chiral phase transitions and an 
assessment of the potential phenomenological observables both in early universe cosmology 
as well as in relativistic heavy ion collisions. We begin this program in this article by focusing 
on the relaxation rate of long-wavelength fluctuations of the order parameter at and near 
the critical point in a consistent non-perturbative framework. 

Strategy: We begin our study of critical slowing down by analyzing the relaxation 
rate of long-wavelength fluctuations of the order parameter at and near the critical point 
in an 0{N) scalar field theory, which is a phenomenological arena to study the relaxation 
of sigma mesons and pions. Our first step is to obtain the relaxation rate to lowest order 
in perturbation theory (two loops). This calculation reveals clearly the breakdown of the 
perturbative expansion for long-wavelength fluctuations at or near the critical point as a 
result of the strong infrared behavior for soft loop momentum and the necessity for a non- 
perturbative treatment. We then implement a non-perturbative resummation of bubble- type 
diagrams via the large N approximation to obtain the damping rate in the next-to-leading 
order in the large limit. The resummation implied by the large limit to order is akin 
to that obtained via the renormalization group with the one loop beta function and reveals 
the softening of the scattering amplitude and the crossover to an effective three dimensional 
theory for momenta q <S AT with A the quartic coupling. 

Summary of main results: 

We have obtained the relaxation rate for long-wavelength fluctuations of the order pa- 
rameter at the critical point and for homogeneous fluctuations near criticality both to lowest 
order in perturbation theory (two loops) and near the critical point to next to leading order 
in the large N limit. 

The two-loop results for the relaxation rate for a fluctuation of wavevector k of the order 
parameter at the critical point is found to be r(/c, T) oc X^T'^/k whereas near the critical 
point, homogeneous fluctuations (with A: = 0) relax with a rate ro(m2-,T) oc X^T'^/rriT- 
Here mr oc \T — Tc\^/'^ <C Tc is the effective thermal mass. These results clearly reveal the 
breakdown of the perturbative expansion in the long wavelength limit A; at T = Tc and 
for T ^ Te and A; = 0. 

A detailed analysis of the different contributions to these results for the relaxation rate 
shows that the rate is dominated by very soft loop momentum q <^ XT which in the weak 
coupling limit A <^ 1 are classical. The implementation of a non-perturbative resummation 
via the large N limit explicitly leads to an effective scattering amplitude that vanishes in the 
static long- wavelength limit as a consequence of the crossover to a three dimensional theory 
for loop momenta q <^ XT. This effective scattering amplitude allows us to recognize that 
the effective three dimensional coupling for soft momenta approaches the three dimensional 
non-trivial (Wilson-Fisher) fixed point in the long- wavelength limit near the critical point. 
The large N resummation for the relaxation rate incorporates this effective three dimensional 
coupling in the spectral density that determines the imaginary part of the retarded self-energy 
for the order parameter. Since the effective three dimensional coupling is driven to its fixed 
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point at long wavelength, the contribution from very soft loop momenta q <^ XT which give 
the strongest infrared behavior in lowest order in perturbation theory is effectively screened 
by this renormalization of the coupling. Consequently the most important contribution to the 
relaxation rate arises both from the semisoft classical region of loop momentum T ^ XT 
and also from the hard region q > T. A detailed analysis of the contribution from the loop 
momenta reveals a non-perturbative ultrasoft scale 
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We find that for soft momenta k ^ k^s the damping rate is dominated by classical semisoft 
loop momenta and given by 
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For k , rriT <C k^s the classical approximation breaks down and the damping rate at the 
critical point mj' = for k <C k^g is given by 



Tik,T) 
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(1.2) 



For homogeneous fluctuations near the critical point {k = 0, rriT oc |T — Td^^^ ^ 0) the 
damping rate is given by 



To{mT,T) 



AnT 



l + O 



In^ 

my 



;i.3) 



Thus critical slowing down, i.e, the vanishing of the quasiparticle width F for long- 
wavelengths emerges in the ultrasoft limit k <^ kus or very near the critical point itlt <^ kus 
where it vanishes logarithmically slow in the k,mT — > limit to this order in 1/A^. Notice 
that in such regimes the rate is independent of the coupling A. 

The large approximation is not limited to weak coupling and our results apply just 
as well to a strong coupling case A > 1 wherein we find that the relaxation rate is given 
by eq. (|1.2|) - (|1.3|) . However this analysis clearly reveals that for weak coupling there emerges 
a hierarchy of widely separated scales for loop momenta: from hard g > T to semisoft 
T ^ q ^ XT, and soft AT ^ q that lead to different contributions to the relaxation 
rate. Which is the relevant scale for the damping rate is determined by the wavevector of 
the fluctuation of the order parameter and the proximity to the critical temperature. For 
k , rriT ^ kus the classical approximation does apply and the damping rate is dominated by 
the soft and semisoft classical loop momenta [with the result ( |1.1| )], whereas for k , my <^ kus 
the classical approximation breaks down and the damping rate is dominated by hard loop 
momenta q >T [with the results (|1.2D-(PT5D]. 
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A similar hierarchy exists in non-abehan plasmas |T8|- |^ and we compare and contrast 
our results in the scalar theory with those in the hard thermal loop approximation in abelian 
and non- abelian plasmas |]18|- . 

This article is organized as follows: in section II we introduce the model, obtain the real- 
time equation of motion for the order parameter and describe the strategy followed to obtain 
the relaxation rate. In section III we carry out a perturbative analysis of the relaxation 
rate to two loops order, recognize the breakdown of perturbation theory and compare to 
the case of the hard thermal loop resummation program in gauge theories. In sections IV 
and V we introduce the large N limit, obtain the effective static scattering amplitude in 
leading order in the large N and discuss the dimensional crossover for soft momenta and 
the effective three dimensional coupling being driven to the three dimensional fixed point. 
We then use these results to obtain the relaxation rate and near criticality to order 0{1/N) 
in the large N limit and explicitly discuss the screening of the soft loop momenta. The 
contribution from classical soft and semisoft momenta and that of hard loop momenta are 
analyzed separately to highlight the important differences. In this section we discuss further 
the validity of a quasiparticle interpretation of the collective long-wavelength fluctuations of 
the order parameter. In section VI we summarize our conclusions and results and discuss 
the next step of the program. In appendix A the equations of motion in the large N limit 
are derived and in appendix B the polarization integral is computed. 



II. PRELIMINARIES: THE MODEL AND THE STRATEGY 

We study the model of scalar fields in the vector representation of 0{N), which is 
conjectured to describe the equilibrium universality class for the chiral phase transition with 
two light quarks for = 4 H]. The Lagrangian density is given by 



2 



A 



C= -{d^^y - - m'T + 6m'{T) - — [<l>^(x)]' + J ■ $ (2.1) 



2A^' 



where the external current J has been introduced to generate an expectation value for the 
scalar field (i.e. the order parameter) by choosing it to be nonzero along a particular (sigma) 
direction. 

The counterterm 5m? {T) is introduced to cancel the tadpole contributions (local terms) 
so that perturbation theory (or the large A^ expansion) is carried out in terms of the effective 
thermal mass mx- In particular to leading order in the large A^ expansion there is the hard 
thermal loop contribution given by the usual tadpole term [^] oc A($^)/A^ oc A which 
combined with the zero temperature (negative) mass squared leads to an effective finite 
temperature mass oc (T^ — T"^). The critical theory corresponds to T = Tc, i.e. mx = 0. 
In this case the counterterm 5m? {T) is adjusted consistently order by order to set the effective 
finite temperature mass equal to zero. 

As stated in the introduction, our goal is to obtain the relaxation rate (damping rate) of 
the order parameter at and near the critical point. This will be achieved by obtaining the 
equation of motion for the expectation value of the scalar field, i.e, the order parameter and 
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treating its evolution in real time as an initial value problem. This is achieved by coupling 
an external source that serves the purpose of preparing the initial state. From the equation 
of motion we recognize the self-energy and compute the relaxation rate from its imaginary 
part on shell. We write 

, t) = ^{x, t) 5^^^ + r]\x, t) ■ (r/(f , t)) = (2.2) 

where we chose the particular direction "1" by choosing external source term to be differ- 
ent from zero along this direction to give the field an expectation value (see below). The 
equation of motion for (p{x,t) is obtained by imposing that {r]^{x,t)) = consistently in the 
perturbative expansion |2^ . In terms of the spatial Fourier transform of the order parameter 



ip and following the steps detailed in appendix A. 2 (see also |Q) we find 

(Pkit) + [k^ + ml + 5m\T) + mL(T)] ^kif) + / ^retAt - t') Mt') dt' = Jk{t) 



where Jk{t) is the external source that generates the initial value problem and Tjret,k{t — t') 
is the two-loops retarded self-energy without the tadpole contributions. The one and two- 
loops tadpole contributions (local) are accounted for in mf^^(T). As described above, the 
counterterm 5m? {T) is fixed consistently in perturbation theory by requesting that it cancels 
all constant (in space and time) contributions to the self-energy (such as the tadpoles) i.e, 

W(T,)+mL(Te)=0 

The retarded self-energy has a dispersive representation in terms of the spectral density 
p{u, k) given by 

- t') = / ^ e-"'<'-''> /d. -^^^ (2.3) 

J Ztt J 00 — 00 — te 



in terms of which the relaxation (damping) rate is given by |18 



= "2^ ^^"^'^^^^ ^' ^^-^^ 

where O0p{k) is the position of the pole in the propagator, i.e, the true dispersion relation. 
For the perturbative two loops or to leading order in the large N limit as studied here 
ujp{k) = k"^ + which at T = Tc takes the form O0p{k) = \k\. 

With the purpose of clearly revealing the breakdown of the perturbative expansion for soft 
momenta k -C AT we will begin our analysis by focusing first on the perturbative evaluation 
of the damping rate. At one loop order the only contribution to the self energy is given by 
the tadpole term A($^(x, t))/A^ which is local, determines to lowest order the temperature 
dependent mass rriT oc |T — T^^/^ and determines the counterterm E^. Furthermore this is 



the leading contribution in the hard thermal loop limit |T8| , |22[] . The lowest order contribution 
to the absorptive (imaginary) part of the self-energy arises at two loops and is studied in 
detail in the next section. 
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III. PERTURBATION THEORY: TWO LOOPS 



We begin our study by carrying out a perturbative evaluation of the damping rate to 
lowest order, i.e. to two loops to reveal several important features of the soft momentum 
limit, and to pave the way to implement a non-perturbative evaluation of the self-energy in 
the large N limit. Furthermore, as it will become clear during the course of the calculation, 
the lowest order contribution contains some of the important ingredients of the large 
limit and will highlight the contribution to the relaxation rate from different regions of loop 
momentum. 

After substracting the one and two loops tadpole contributions which are cancelled by 
the counterterm, the spatial Fourier transform of the retarded self-energy reads: 

E„,,,(i - (') = 8A^ ^ / {Gljt - f) G^t - (') G|(i - *') 

- Gijt - t') Gl^,4.t - f) G}{t - f)} e(t - f) (3.1) 

where the Wightmann functions G^,G^ are given in Appendix A.l. 

With the purpose of comparing with the results of later sections, it proves convenient to 
introduce the intermediate quantities 

- t') = -2A ^ / ^ GUt - 1') G}{t - t') 

= j dqoe-'^'^'-''^ S>{qo,q) 

- = -|(iV + 2) / ^ GUt - t') G^{t - t') 

and using the expression for the Wightmann functions G^if.'-t'). G^{t — t') given in appendix 
A.l, it is a straightforward exercise to show that the spectral functions S^{qo,q); S^{qo,q) 
obey the KMS condition 

5<(?o,g) = e-^«°5>(go,g) 

Introducing the spectral density 

(^{Qo, q) = S>{qo, q) - S<{qo, q) 
which at the critical point T — Tc, i.e. rriT = is found to be given by 

jV + 2 r (Pp 1 

a{qo, q) - 2A — ^ j J^nf^p^pT^^^^ ^ '^'^^^^ '^^^ ^'^^^^ - |P + ?1 - P) - <^(?o +\p + q\+p)] 

+ [np - n^+p] - |P + ?1 + P) - K% + |P + ?1 - P)]} (3.2) 

we find 



9 



S^iqo, g) = [1 + n{qo)] a{qo, q) , 5'<(go, q) = n{qo) a{qo, q) 



(3.3) 



where 



(3.4) 



The near critical case is considered in sec. III.C. 

A lengthy but straighforward calculation with the Bose-Einstein distribution functions 
for massless particles leads to the following expression for a (go, q) 



c^(go, g) = g 2 < ©( go - g) sign(go) + — In 
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(3.5) 



It is important to emphasize that the second, finite temperature term is the combination of 
two different contributions: a two massless particle cut with support in the region go > g 
and a Landau damping cut with support in the region — g < go < g. 
Introducing the Fourier representation of the theta function 



e(t-t') = -/ 



27rz cij + ze 



(3.6) 



in the expression for the self-energy ( |3.1| ), we find the spectral density that enters in the 
dispersive representation of the retarded self-energy (|2.3|) to be given by 



piuj,k) 



4A f d^q dqo 



N J (27r)3 2|fc + g-1 "^--U^ ■ 



o-(go,g) 1 + «o-+fc + '^(go) 5(tu-go- |/c + g1) -5(c^ + go + |/^ + gl) 



AX r d^q 1 r 



- |A? + g1) [n.^g - n(|A^ + g] - ^) 



iV^ (27r)3 2|^+g-1 
+ o-(g,t^ + |A; + gl) - n(u; + |A; + g])]} 

with o"(go, g) given by ( p3|) and we used that, 

cr(-go, g) = -o-(go, g) , i + ^(-go) = -?^(go) • 

Although this seems to be a cumbersome manner to write down the two loop contribution, it 
does prove convenient to establish contact with the large N description in the next section. 

We are interested in the relaxation of long- wavelength fluctuations of the order parameter, 
hence we will consider the soft limit k <^T. 

We study the contributions coming from soft q <^ T and hard q > T loop-momenta 
separately. 



(3.7) 



10 



A. The soft momenta contribution (classical region) 



This is the classical region where the Bose-Einstein distribution functions can be approx- 
imated as n{uj) ^ T/uj, ?a T /k. 

In this regime the contribution of the soft momenta q <^T yields 



27r2 



iV2 



In 
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\k + q\ 



+ 



LLi+\k+q\+q 
Ld+\k+q\—q 



\k + q\ 



where we only kept the contribution of order T to cr(go, q)- We evaluate the spectral density 
aX uj = k leading to the following form for the soft momenta contribution to the damping 
rate (O) 



Td{k,T) 




N + 2 



[Mk,T) + Mk,T)] 



with Ji(fc, T) and J2{k,T) given by the following expressions 



(3.8) 



Mk,T)= d\ 



J2{k,T)= d'q 



In 
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q\k + q\'^ {\k + q\ + k) 



(3.9) 



The angular integrals can be performed analytically and k can be scaled out of the integral 
by introducing the variable x = q/k leading to the final expression 



Tci{k,T) = 
with the function F[x\ given by 

F[x\ = 



X^T^ N + 2 
IQn^k m 



dx F[x\ 



(3.10) 
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and is depicted in Fig. 1. In principle, the upper limit in the integral in (|3.10| ) should be aT / k 
with a -C 1 to restrict the integral to the soft momenta where the classical approximation 
is valid, but the integrand falls off as 1/x^ for x ^ 1 and the integral is dominated by the 
small X region < x < 1 as shown in fig. 1. The integral from x = to x = oo in eq.( ^.l(J| ) 
gives 27r^ and the contribution from the classical loop momenta q <^T to the damping rate 
to two loops order is thus given by 



X^T^N + 2 
'8^ W 



(3.11) 
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As it will be seen in detail in the next section, the large N resummation leads to a screening 
of the soft loop momentum which cuts off the contribution of momentum q < XT. Hence, 
with the purpose of comparing the perturbative two-loop result with that in the large N 
limit, it proves convenient to obtain the contribution to the damping rate from the region 
of semisoft loop momenta with AT < q < aT with A ^ a ^ 1 for the case of soft external 
momentum XT/k ^ 1. 

The contribution to the integral of the function F[x] from this region is given by 



therefore this region of loop momentum gives a contribution to the damping rate given by 

AT AT + 2 AT 



T>g>AT 



27T m k 



In—. (3.12) 



Thus we see that the contribution from the soft region of internal loop momentum q < XT 
contributes a factor XT/k ^ 1 larger than the region of momentum q > XT for soft external 
momentum k <^ AT. This observation will become important when we compare with the 
large result, because as we will show explicitly below, the resummation of the effective 
scattering amplitude will lead to a softening of the effective vertex and hence screening of 
very soft momenta in the loop. 



B. The hard momenta contribution 

We now focus on obtaining the contribution to the damping rate from hard loop momenta 
q > T. From the expression of the finite temperature contribution to the spectral density 
cr(go, q) ( |3.5| ) it is clear that hard momenta qo,q >T will be exponentially suppressed unless 
either |g — go| -C T or |g + q'o| -C T. Consider the expression for the spectral density (IsTTI) for 
uj = k and consider the contribution from the delta function with support for go = k — \k + q\\ 
it is straightforward to see that the other delta function will give a similar contribution. For 
q > T ^ k we find that |go + qI = k\l ~ cos6\ where 6 is the angle between k and q 
and \qo — q\ > 2T, hence the region of loop momentum that dominates corresponds to the 
emission (or absorption) of a pair of scalars (the particles in the loop) with total center of 
mass momentum collinear with the external momentum. 

Keeping the leading term oc ln(|go + '?|/2T) in (T(go, <?) and the full occupation factors in 
the expression for p the spectral density ( |3.7] ) becomes 

T 

pioj = k,k) (X — A^ T In — 

k 

and the contribution to the damping rate from the hard loop momentum region is 

T 

Thardik,T) oc A^ T In — 
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which is a factor (k/T)\nj ^ 1 smaller than the contribution from the classical loop 
momenta (|3.11|) for k/T <^ 1. In summary, when the external momentum is soft k <^ XT 
the damping rate is completely determined by the classical region g <C T of loop momenta 
and given byeq. fpID . That in a scalar theory the leading temperature effects are determined 
by the classical region of loop momenta was already anticipated in references p^ , |25| but the 
computation above identifies the contribution to the relaxation rate from several different 
regions of loop momentum. This identification will become important to understand the 
result obtained from the large limit. 



C. Near criticality 

A calculation very similar to that in the critical inhomogeneous case can be carried out 
for homogeneous fluctuations {k = 0) near the critical point with the effective thermal mass 
m"^ oc A(T^ — T^) straightforwardly. In this case the angular integrals are trivial and most 



of the steps are similar to the critical case leading to (see also p4 

X'^T^ N + 2 



To{mT,T) 



A similar analysis for contributions from different regions of loop momentum is obtained by 
replacing k rriT in the arguments above. 

The resonance parameter T{k,T) / ujp with ujp the position of the single (quasi) particle 
pole determines how broad is the resonance. If r(/c, T) <^ ojp the quasiparticle can be 
described by a narrow resonance and its decay occurs on time scales much longer than those 
of the microscopic oscillations On the other hand for T{k,T) ^ ujp the notion of 

quasiparticle is not appropriate and the excitation is described by a very short lived broad 
resonance. 

The two loops calculation reveals that at the critical point uit = 

T{k,T) + 2 (XtV 

k SnN^ \ k J ' 

analogously, near criticality for homogeneous fluctuations k = 

T{mT,T) T>^^ N + 2 /Xr_y 
rriT 8 ttA^^ \mT J 

Hence up to this order a quasiparticle interpretation is not reliable for k , rriT <^ AT. 

Moreover, for very soft external momenta or very near the critical temperature, k ; -C 
AT the perturbative expansion clearly breaks down and a non-perturbative scheme must be 
invoked to study the damping rate. 

This situation is similar to that in the hard thermal loops (HTL) program in the sense 
that for external momenta k -C AT (in gauge theories A must be replaced by the gauge 
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coupling squared) a non-perturbative resummation is needed [0,^. However, here the 
similarity ends and the major difference with the HTL program is revealed: whereas in the 
HTL case the non-perturbative region is dominated by hard internal loop momenta q >T, 
the relaxational dynamics at the critical point is dominated by soft classical (and as it will 
become clear below, also semi-soft) internal loop momenta q <^T . The difference can also 
be clearly seen formally by restoring the h in the contributions: the temperature always 
appears in the combination T /h (from the distribution functions), in the HTL program the 
gauge coupling constant squared — >■ e^h (since this is the loop counting parameter) hence 
the HTL scale e^T^ — > e'^T'^/h. However in the scalar case the loop counting parameter is 
A Xh, hence the contribution A^T^ is classical i.e. independent of h. Therefore whereas 
in the HTL program perturbation theory breaks down at a semiclassical scale k oc eT/^/h, 
at the critical point of a scalar field theory the perturbative expansion breaks down at a 
classical scale k oc AT. In the HTL program the damping rate of collective excitations 
is typically of order e^T and the quasiparticle poles (plasmons and plasminos) are of order 
Up oc eT/ \/h hence for weak coupling the long-wavelength quasiparticles are always relatively 
narrow resonances. This is in striking contrast with the case of a critical scalar theory where 
the long-wavelength excitation of the order parameter is gapless. 



IV. LARGE N 

Having recognized the non-perturbative nature of the relaxation for the long-wavelength 
components of the order parameter, we seek to use a consistent non-perturbative description 
and study the relaxation of the order parameter in the large limit. This limit is best 
studied by introducing an auxiliary field that replaces the quartic interaction via a Gaus- 
sian integration (Hubbard-Stratonovich transformation), hence the Lagrangian density 
becomes 

^^{x) + ^X^{x) + J{x) -^ix) . (4.1) 

Before we engage in a study of the damping rate, it is important to highlight that the 
large N expansion effectively provides a reorganization of the perturbative series which 
for example at leading order and at zero temperature is akin to the resummation of the 
leading logarithms via the renormalization group for the scattering amplitude. We now 
study in detail this resummation at finite temperature which will reveal the screening of the 
scattering amplitude for soft momenta which in turn will be responsible for screening the 
infrared behavior of the spectral functions and the damping rate. 

In this section we will focus on the critical theory rriT = 0. The analysis of the off-critical 
case is given in section V. 
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A. Effective Scattering Amplitude 



To leading order in the large the two particle to two particle scattering amplitude is 
dominated by s-channel exchange and is completely determined by the propagator of the 
auxiliary field x- Fig- 2 shows the Dyson sum for the propagator of the auxiliary field in 
leading order in the large limit and fig. 3 shows the s-channel scattering amplitude in 
leading order, the t and w-channel contributions are subleading. The bubble diagram which 
is the building block of the propagator of the auxiliary field (and therefore the s-channel 
scattering amplitude) is simpler to be calculated in the Matsubara formulation of finite 
temperature field theory with an external frequency Vn = '^irnT and given by 

/,.K,.)^2ATE/|^ ^J-.,„^^^^;^,.^^-,. . (4.2) 

To illustrate the resummation in a more clear manner, we focus on the static limit which 
is obtained by setting the external Matsubara frequency to zero. The strongest infrared 
behavior and leading contribution in the high temperature limit arises from the term m = 
in the Matsubara sum, the remaining spatial momentum integral is carried out leading to 

AT 

hub{0,q) = — (4.3) 
Aq 

and the s-channel scattering amplitude in the static limit is given by (see Fig. 3) 



2' 2' iVl + /,„,(0,g) 

We see that the effective temperature and momentum dependent coupling constant, defined 
as the coefficient of in the s-channel scattering amplitude in the static limit, is given by 

^e//(g) = • (4.4) 

This expression reveals several noteworthy features. Firstly, we see that at high temperature 
in the critical region the actual expansion parameter in the sum of bubbles is XT/q with q 
the spatial momentum tranferred into the loop. The factor T is a consequence of the dimen- 
sional reduction and the factor AT can be interpreted as the dimensionful three dimensional 
coupling. Since the expansion is in terms of dimensionless quantities the factor q in the de- 
nominator is required for dimensional reasons. In fact this can be understood via a parallel 
with the calculation at zero temperature in 4 — e space-time euclidean dimensions with a 
coupling AT*^ with T now some dimensionful scale, the loop integral for the massless theory 
produces a factor and for e = 1 i.e. the three dimensional theory one finds the result 
for the finite temperature loop in the static limit. Secondly, the expression for the effective 
coupling ( [4. 41) is a result of the large N resummation to leading order and is the same as that 
obtained from the solution to the renormalization group equation for the running coupling 
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using the one-loop beta function obtained in the e expansion and setting e = 1, i.e. the 
large N resummation is akin to the resummation obtained from the renormalization group 
in euclidean field theory, in the sense that the leading order in the large leads to a run- 
ning coupling which is the same as that obtained from the one-loop beta function. Thirdly, 
since the effective expansion parameter in the sum of bubbles is XT/4q it is convenient to 
introduce the three dimensional coupling A3(g) = XT/iq and its effective counterpart 



Xsjq) _ AT 
1 + A3(g) Aq + XT 



Keff{q) = = ■ (4-5) 



The main point is that this effective three-dimensional coupling is driven to the three- 
dimensional Wilson-Fisher fixed point A* = 1 in the soft momentum limit g — > 0, while 
the effective four dimensional coupling (|4.4|) is driven to the trivial fixed point in this limit. 
Hence, whereas the three dimensional coupling A3(g) = XT/Aq diverges in the q ^ limit, 
the large N (equivalent to the renormalization group) effective coupling X^^effil) is driven 
to a finite fixed point in the soft momentum limit. Therefore the large resummation is 
effectively screening the infrared divergences associated with the soft momentum limit much 
in the same manner as the resummation implied by the renormalization group within the 
e expansion. Obviously, in exactly three euclidean dimensions one could hardly justify the 
validity of an e expansion, but the large limit provides a non-perturbative framework that 
includes a similar resummation. The main point of this discussion is the realization that 
the resummation implied by the large N limit provides an effective coupling constant that is 
well behaved in the infrared limit, thus leading to the conclusion that the simple point-like 
scattering vertex must be resummed before attempting to compute the damping rate or any 
other transport coefficient near the critical region. 

The analysis in this section reveals the role played by the scale AT: internal loop momenta 
q <^ AT lead to non-perturbative contributions, in the weak coupling limit A <^ 1 these non- 
perturbative scales are classical, on the other hand for q ^ AT the effective couplings (either 
four or three dimensional) are small for weak coupling A and the effective vertices coincide 
with the bare vertices. The implications of this discussion will be important to understand 
the different contributions to the relaxation rate. 



B. The relaxation rate 

As discussed above at leading order in the large A^ limit the only contribution to the 
scalar self-energy is a tadpole oc A($^)/A^ ~ C^(l)) which results in the effective thermal 
mass ruT oc |T — Td^/^ and is cancelled by the mass counterterm. In this section we consider 
the theory at the critical temperature where the renormalized temperature dependent mass 
exactly vanishes. 

At next-to-leading order 0{1/N) the self-energy obtains an absorptive part and is given 
by the diagram shown in fig. 4. 

In appendix A. 2 we provide the details necessary to obtain the retarded self-energy in 
terms of a dispersion relation as in eq.( p.3| ), with the spectral density 
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4A 



NJ (27r)3 2|A; + g-1 



2A 



P(?o,g) + n^+p + n{qo)] 6{uj - qo - \k + ^) - 6{u + Qo + \k + q\] 



(Pq 



[p{uj- \k + q\,q) 



n 



NJ (27r)3|^ + gj 
+ P(^ + \k + q],q) [n^_^^-n{uj+ \k + q])]} 



q+k 



- n{\k + q\- oj) 



We performed here the integral over go by using the delta functions thereby setting the 
combination qo + \q + k\ = ±a; for the respective delta functions. 
In appendix A.l we show in detail that 



p{(lo,(l) 



n/(go,g) 



TT [l + URiqo,q)Y + njiqo,q) 



(4.7) 



where Ili{qQ,q) is given by the leading order in the large N limit of the two loop spectral 
density (|3.2|) , as 



n,(go,g) = 2A7r I 



1 



(27r)3 4p|p + g1 



{[1 + ng-+p-+ rip-] [5(go - \p + q\ - p) - 5{qo + \p + q\ + p)] 
+ [up - n^+p] [5(go - |p + gl + p) - 5(go + |p + gl - P)]] ■ 



The first term, proportional to the sum of the occupation factors, corresponds to the two 
particle cut while the second term proportional to the difference is obviously only present in 
the medium and corresponds to Landau damping [l^. The real and imaginary parts of the 
polarization of the auxiliary field n(q'O) are related by a dispersion relation, i.e. 



If 1 
ni?(go, q) = - duj Ui{q, u) V 

TT J UO — 



(4.8) 



Keeping the leading temperature dependence we obtain 

1 



AT 

n/(go,g) = — In 

Airq 



\io+i\ 
e 2T 

. ko-g| 

i — e 2T 



+ o(at°) 



(4.9) 



In appendix B we show explicitly that the leading temperature dependence for the real part 
of the polarization operator of the auxiliary field is given by 



AT 

nR(go, q) = ^ [e(g - go) - e(-g - go)] + O {X InT) 



(4.10) 



which in the static limit reduces to ( [4.31 ). 

It is clear from eq. (|4.9|) that just like in the case of perturbation theory up to two loops, 
there are two important regions to consider: i) the classical region with go, g <^ T and ii) the 
hard region with go, g > T but with either |g — go| ^ T or |g + go| <^ T the other regions of 



(4.6) 
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hard momentum being exponentially suppressed. We will be primarily interested in the case 
of soft external momentum k -C AT i.e. long-wavelength fluctuations of the order parameter. 
In the high temperature limit the spectral density p{qo, q) takes the explicit form 



p{k+ \k + ^,q) 



p{k- \k + q\,q) 



Aq T 



i2. + 1 
AT ^ 



AT 



^ 2 



+ L2 



where we used eqs.([4.9|)-( |4.10| ) analyzing carefully the support of the theta functions in 
n_R(go) q) and defined 

k±\k + q\+ q 



L± = log 

k ± \k + q\ — q 

Introducing in the integral ( ^.61 ) the dimensionless variable y = q/k and x = cos6 (where 6 
is the angle between the vectors k and q) and setting uj = k yields for the damping rate (|2.4| ) 



T{k,T) 



fc2 



-K N T Jo 



y^dy 



+1 dx 



+ 



1 



1 w{x,y) 
1 



^Pkw(x,y) _ I ^0k[w{x,y)-l] _ ]^ 

L-{x,y) 



L+{x,y) 



AT 



where 



oP k w{x,y) 



L±ix,y) = log 



of3k[w{x,y)+l] 



lTw{x,y) + y 



AT 



+ Ll{x,y) 



Ll{x,y) 



(4.11) 



and w{x, y) = J 1 + y'^ + 2xy . 



lTw{x,y) -y 

For XT/Ak ^ 1 the region of small y (small momentum) is screened by the resummation of 
the scattering amplitude and leads to a small contribution to the damping rate of order k. 
For y > XT/k the screening is not effective and the integrals in eq. (|4.11|) are dominated by 
the neighborhood of the point y = y* aX which L\{x,y) is of the same order as the other 
square in the denominators, i.e, [4/cy/AT]^. For large y we can expand w{x,y) as follows 



w{x, y) = y + X + O 



(4.12) 



We thus find 



y 



AT , AT 
m 



4:71 k Atx k 

[For k ^ T the 1 in the denominator of the first term of eq. ([4.11| ) is irrelevant in the 
determination of y*] . 



18 



The classical approximation for the occupation numbers, i.e, loop momenta <^ T can be 
used at y = y* provided k y* <^T and A ^ 1, i.e, for k ^ kus where 

A T _47r 

kus = -;— e ^ 

is the ultrasoft scale. 

As it will be discussed in detail below, for k < k^s hard momenta dominate the contri- 
butions to the width r(fc, T). 

We analyze in the subsequent section the width T{k,T) in the two regimes : soft for 
which T ^ k ^ kus and ultrasoft k <C kus- 



C. Classical contribution 



To obtain the contribution from classical momenta q <^ T we perform the following 
approximations: i) approximate 11/(^0, q) by its limit for qo,q <^T 



AT 



Qo-Q 



with the real part given by the leading temperature contribution ([4.10|) and ii) approximate 
the Bose-Einstein occupation factors by their classical limit, i.e. in (|4.6|) we replace 



1 + ng-+p-+ n(go) = T 



Qo\q + p\ 

We thus find from eg. ( ^.111 ) for the classical limit contribution to the damping rate. 

k ^T/k „ /•+! rfx I 1 L+(x,y) 
y dy ' 



+ 



7lN Jo 

1 



1 w'^{x,y) \ l-w{x,y) ^^-j-^ + Ll{x,y) 



L-{x,y) 



1 + w{x, y) 



+ LUx,y) 



(4.13) 



with X^^effio) being the effective three dimensional coupling given by ( |4.5| ) and A* = 1 
the three dimensional fixed point. We have introduced an explicit upper momentum cutoff 
(Imax = (y-T with a <^ 1 that restricts the integration domain to the region where the classical 
approximation is valid. 

The expression ( [4.13[ ) clearly reveals the role of the three dimensional effective coupling 
X^ effiq) given by ( |4.5| ) and its non-trivial (three dimensional Wilson-Fisher) fixed point 
A* = 1 reached in the soft limit g — > 0. The phenomenon of screening of the infrared 
behavior by the renormalization of the coupling is now explicit, the region of soft loop 
momentum q ^ AT is independent of the coupling and temperature because the effective 
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three dimensional coupling is near its non-trivial fixed point. The only scale in the integral 
in the soft-momentum region is k and a dimensional analysis reveals that the contribution 
from this region is proportional to k. 

On the other hand, when the loop momentum is g ^ AT the renormalization of the cou- 
pling is ineffective and the effective coupling coincides with the three dimensional coupling 
XT/Aq. For weak coupling A <^ 1 and loop momenta q ^ AT the effective three dimensional 
coupling is Xz^effil) ~ A3(g) = XT /Aq <^ 1 and the denominators in eq.( [4.13|) are dominated 
by the terms 1/ X^^effio)- If the logarithms can be neglected we clearly see that this contri- 



bution is the same as that given by the integrals Ji(fc, T) and J2(^, T) given by eqs. (|3.9|) in 
the two loop computation of the damping rate ( |3.8| ), which is proportional to X^T"^ /k. This 
region begins to dominate for q > XTln{XT/k) when l/X^^effi^) becomes larger than the 
logarithm in the denominators in both integrals in eq.( |4.13| ). If the loop momenta are such 



that T ^ g in this region, i.e. the classical approximation is valid, the contribution of this 
region to the integral can be estimated by cutting off the integrals in eq. ( [4. 131 ) at a lower 



momentum of order g^m ~ AT ln(AT/ k). Hence following the same arguments as for the two 
loops case that led to the estimate we conclude that the region of classical semi-soft 

loop momentum T ^ q^ qmm leads to a contribution to the damping rate ~ AT. However 
as k becomes smaller, g^m approaches the cutoff aT, i.e, the limit of validity of the classical 
approximation and the logarithmic terms cannot be neglected. In particular for k < kus 
with kus the ultrasoft scale introduced above the integral becomes sensitive to momenta of 
order of T and the classical approximation breaks down. For these ultrasoft momenta of the 
fluctuations the damping rate is determined by the region of hard loop momentum g > T. 

This analysis yields to a preliminary assessment of how different regions of loop momen- 
tum will contribute to the damping rate: 

i) the soft region of loop momentum g <^ AT is dominated by the three dimensional fixed 
point and contributes to the damping rate 



g<AT 



OC k 



ii) the semisoft region of loop momentum T ^ q ^ XT is still dominated by classical modes 
but the renormalization of the scattering amplitude is irrelevant. If qmm ~ ATln(AT/A;) <C 
T the logarithms can be neglected and the integrals in eq.( [4.13D behave similarly to the 



perturbative two loops computation. For momenta T ^ g ^ qmin ~ ATln(AT/A;) the 
integrals are dominated by the terms in the denominators proportional to 1/A|(g) leading 
to a contribution 

rT><7>g„i„ OC AT 

The validity of the classical approximation and the dominance of semisoft loop momenta is 
warranted for weak coupling when there is a clear separation between the hard scales with 
g > T the semisoft scales with T ^ g > AT and the soft scales for which AT ^ g. 

iii) However when q^m > T the logarithmic terms are the dominant terms in the de- 
nominators of the integrands and this region is sensitive to the hard loop momenta g > T 
and the classical approximation is not warranted. This region requires the full Bose-Einstein 
distributions and will be studied in detail below. 
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After this preliminary assessment, we now provide a quantitative analysis of the different 
regions. 

As argued above in the region y <^ A3(A;) = \T/{4:k), i.e. q -C AT, the integrals in 
eq. (^4.13|) are independent of A3 (A;) (i.e, of both the coupling constant and temperature) and 
infrared finite, contributing to the damping rate a term that is proportional to k. This is the 
region of loop momenta for which the effective coupling is near the three dimensional Wilson- 
Fisher fixed point. The screening of the scattering amplitude is ineffective for loop momenta 
q > XT i.e, the semisoft scales, in this region y ^ Xsik) lnA3(fc) and the term [yX^i^k)]'^ 
dominates over the logarithms leading to a contribution to the damping rate proportional to 
kX^ik) = XT/4:. Since in this region y > Xs{k) ^ 1 we can approximate w{x,y) according 
to eq. (|4.12| ) and the integrals simplify considerably. 

Up to l/X-i{k) corrections we find 



rd(fc,T) 



k f^T/k 



+ 



Nn Jx3(k) 
In^ 

l+x 



dy 



+1 



dx 



In^ 

l—x 



y 



In^ 



2 2u 



1 + 



A3(fe) 



+ 1 



+ ln2^ 



(4.14) 



Now the angular integrals (over the variable x) can be performed changing the integration 
variables to y = A3(fc)u and expanding in inverse powers of lnA3(fc), since for A ^ a ^ 1 



In 



2A3(A;)m' 
l±x 



ln[2A3(A;)] 



1 + 



MHk), 



This is certainly a slowly converging approximation for soft momenta but numerical calcu- 
lations show it to be reliable (see figs. 5-7 below). Then we find from eq. (|4.14|) 



rci{k,T) 



2k Xsik) 

Nn 
XT 
2Nn 



du 



ln[2A3(fc)] 



+ 



ln[2A3(fc)] 



1 + 



7r2(« + 1)2 + ln^[2A3(A;)] (nu)^ + ln^[2A3(A;)] J 
1 ^ 



1 + 



MHk), 



InX^ik)^ 



A detailed analysis of both integrals above reveal the presence of two important scales: a) 
the cutoff scale y = aT/k, that is g = aT with a -C 1 that determines the regime of validity 
of the classical approximation when the Bose-Einstein occupation factors may be replaced 
by their classical counterparts and b) a scale y* ^ Xs{k) InXslk) + ■ ■ ■ at which there is a 
crossover of behavior in the denominators of the integrals. For y < y*, y/X^{k) <C Iny and 
the denominator is dominated by the logarithm, whereas for y ^ y*, y/X^{k) ^ Iny and 
the integrands behave as Iny/y"^ which is the same behavior as that in the integrals Ji{k, T) 
and J2{k,T) in eqs.( ^.9D for the two loop computation. In the case when y* <^ aT/k the 
integrand falls off very fast and the integral is independent of the upper cutoff. 

The result ( [4.15| ) is confirmed by a careful numerical study of the integrals in this range 
and displayed in figs. 5-7. 



(4.15) 
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The condition that y* <^ aT/k translates into the following condition for k 

k > ATe"^ 

It is clear that y* becomes of the order of T /k and therefore the crossover scale becomes of 
order T for /c ~ kus — ^ . Hence for wavevectors k ^ k^s the classical approximation 

is valid and the damping rate is dominated by semisoft classical loop momenta T ^ q ^ XT 
and given by 



l + O 




XT 4 

for k > e^^ . (4.16) 

An 



In the opposite limit, i.e. for k -C kus the crossover scale y* ^ aT/k and the terms 
[y/Xslk)]'^ in the denominators are negligible as compared to (Iny)^ in this range. In this 
case the integrals can be evaluated by neglecting the [y/X^lk)]'^ in the denominators with 
the result 

This cutoff dependence signals the breakdown of the classical approximation since the inte- 
grand is sensitive to hard momenta of order q > aT. For k -C kus the crossover scale y* 
becomes of order T/k, i.e, q T and we must keep the full occupation numbers, this is the 
regime dominated by the hard loop momentum, which is studied below. 

Thus we conclude that the non-perturbative region of wavevectors for which the classical 
approximation is valid is AT ^ k ^ ^e^^ and in this region the relaxation rate is given 

4n 



by ( [4.16| ). However for long- wavelength fluctuations with wavevectors k -C kus ~ ^ 



the classical approximation breaks down and we must consider the contribution from the 
hard loop momenta. 

This analysis of the classical contribution reveals that i) the screening of loop momenta 
q ^ AT by the infrared renormalization of the scattering amplitude makes the damping rate 
a factor k/XT <C 1 sma//er than the lowest order (two loops) computation, ii) the damping 
rate is independent of momentum for k 3> ^e"^'^/'*' and given by en. ( [4.16|) i.e. there is no 
critical slowing down in the regime of validity of the classical approximation to this order in 
the large N expansion. 



D. Ultrasoft Scale: k ^ e'^'^l^ 



We now focus on the computation of the damping rate in the regime of ultra-soft fluctu- 
ations of the order parameter, i.e, k -C ^ e~^'^l^ . 

In this limit we expand the difference of the occupation numbers to order ^ inside the 
integrand in eq.( ^4.11 ) 

1 \ k e^'='"(^'S') (k\^ 

Tt^ ,., 



Ql3kw(x,y) _ I ^I3k[w{x,y)^l] _ I ~^ J' UPkw{x,y) _ 1)^ \T 
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Since in this region the logarithm dominates, we neglect the terms Anq/XT in the denomi- 
nators. We thus find 



T{k,T) 



NtiT^ 



dy 



'1 dx 



+ 



-1 w{x,y) (^ef^kw{x,y) _ ly 
1 



X 



In order to perform the integration it is convenient to change variables to 



k 

V = -[w{x,y) - 1] 



cr = 



2k w{x, y) — I 
T y + 1- w{x,y) 



The width then takes the form 



r(fc,T) 



T>fc 



+ 



2T 



Nn 



~ e^v^dv da 1^ (l + ^) 



°° e"" dv r°° da 



2k/T (e'" -1) Jv a^ In 



a' ^2 + (i + ^) 
1 



^-1 



We can further approximate these expressions by expanding in inverse powers of Inf. We 



set, 



da ln(l + ^) _ r^da 1 



V + In^ (l + ^) 

da 1 ] 



a^ In^ 



l + O 



In 



T 



V In 



k I 



l + O 



In 



k L 



l + O 



Ml 



The asymptotic form of the width thus becomes 



T{k,T) 



k I 



l + O 



In 



fc / J 



(4.17) 



where we used the integral p? 



°° e" v"^ dv TT^ 

(e- _ ~ y 



There are two important noteworthy features of this result: i) the damping rate for 
ultrasoft fiuctuations is independent of the coupling and ii) critical slowing down of long- 
wavelength fiuctuations emerges in the ultrasoft momentum limit with the damping rate 
vanishing only logarithmically as A; — > 0. 
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The intermediate regime between the soft and ultrasoft scales is difficult to study analyt- 
ically, we therefore studied the damping rate in a wide range of momentum k numerically. 

Figs. 5-7 display the dimensionless ratio NT{k,T)/T as a function of In^ for three 
fixed values of the coupling: A = 12.0, 1.0 and 0.2 respectively as obtained via a numerical 
integration of eq. ( |4.11| ) . We see that the damping rate is a monotonically decreasing function 
of 1^ for large enough values of j. The smaller is the coupling, the slower NT{k,T)/T 
decreases as a function of ^. Furthermore we have established numerically the reliability of 
the results in the soft and ultrasoft regimes, thus conffiming our detailed analysis in these 
cases. 



V. NEAR CRITICAL REGIME: K = 0, \T - Tc\ <C Tc 

Having understood in detail the critical case we are now in position to complete the study 
of relaxation by considering the near-critical case, i.e. |T — -C Tc. The general case of 
k ^ 0,T ^ Tc is rather complicated but we can learn much by focusing on the homogeneous 
case k = 0. There are two important modifications of the previous results that are required 
to study in detail the near-critical case: 

i:) To leading order in the large limit, the finite temperature effective mass squared 
is given by 

m| oc A (T^ - T^) 

therefore the effective mass (inverse of the correlation length) vanishes near the critical 
temperature as m-p oc |T — T^^/^ which is the mean-field behavior consistent. Since the 
absorptive part of the self-energy is next to leading order 0{1/N) we consistently use this 
effective mass near the critical point. Therefore to this order the frequencies are given by 

ii:) The effective static scattering amplitude can be obtained by replacing the massless 
Matsubara propagators in ( [4. 21 ) by the corresponding massive ones, and is now given by 

T \ A 



KrriT rriTJ i + -AHarctE 

2-Kq c 



2mT 



which now reveals the vanishing of the effective coupling 



, , T \ A 

A.ff 



) = (5-1) 



for rriT [compare with eq. ([4.4| 



The critical region of relevance corresponds to m-r <C AT for T — T+. In the case 
under consideration, for homogeneous fluctuations of the order parameter [k = 0) the only 
dimensionful quantity is the effective mass which regulates the infrared behavior of 
the integrals. Thus, just as in the critical case studied above two different regimes emerge 
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which we refer to as: i) the semicritical regime mj^ <^ AT/(47r); ii) the ultracritical regime 
rriT <C AT/(47r)e~^'^''^. It will become clear below that the semicritical and the ultracritical 
regimes correspond respectively to the soft and the ultrasoft regimes discussed at T = Tc, 
k ^ 0. Since the relevant loop momenta are semisoft, T S> g S> AT and hard g > T it is 
clear that the effective coupling ( |5.1| ) behaves just as in the critical case for 7^ studied 
above since for this range of loop momenta q/mx ^ 1. 

In order to compute the damping rate we need the general expression for the resummed 
spectral density p(go; in presence of a non-zero thermal mass, which is now given by [see 
eq.O] 



~( i-\ 4A 



4A 



di^q dqo 



d^q 1 



p{qo, q) [1 + + ^(go)] [K^ - go - ^\k+q\) - 5(w + go + ^\k+q) 



\q+k\ 



\k+q\ 



where ul = 771^ + k'^, n{qo) is the Bose-Einstein distribution function ( |3.4| ) and p{uj±Lj^^^^, q) 



is the massive spectral density at leading order in the expansion, given by equation 
with the following expression for Iii{uj ± "-^1^+^, q) 



Arn^ 



XT , 

sgnfcij) + - — In 




-f3uip 



(5.2) 



We notice that ( p.2[ ) reduces to ( |4.9| ) in the critical limit uit/T 0. 

Keeping the leading correction in m\ yields in the high temperature regime. 



XT 



+0(m^)} 



uj + p 



uo — p 



Arn^ p^ sgn(ci;) 



2^2 



+ 



{uj^ - p" - Ami) - e{uj^ - p") 



log 


uo +p 




uo ~ p 



(5.3) 



The real part ^.^{uo^p) can now be obtained via the dispersion relation 



nij(go, q) 



XT 



Aq 

+0(m^) 



e(g-go)-0(-g-go) + 



Ami 



7t^{q^ - qi) 



mj^ 



.g go. 



(5.4) 
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The damping rate for homogeneous configurations is given by 

7rp(u;fc, k) TTp{mT, 0) 



ro(m^,T) = -lim. 



2mi 



which is now given exphcitly by 

ttA r d?q 1 



ro(mT,T) 



tutN J (27r)3cjq 



{[n{u;q) - n{ujq - m^)] p{mT - ujg, q) 
+ [n{ujq) - n{ujq + uit)] p{mT + ujg, q)} 



(5.5) 



The integral ( p. 51 ) is much simpler than the analogous expression for /c 7^ 0, since the 
angular integration is trivial. Nevertheless, a complete evalution of ( p.5| ) requires a numerical 
integration. We refer to figures 8 and 9 for a numerical evalutation of the dimensionless ratio 
Nro{mT,T)/T in the intermediate (A = 1.0) and small coupling regimes (A = 0.2). 

Just as in the critical case in the near-critical regime with rriT -C AT/(47r) the integral 



(5.5) is dominated by loop momenta of order q > q* 



XT 



XT 



^ ^ ^ rriT hence we 



Ait 



can approximate 117(50, q) and Ilji{qQ, q) at go = ^t ± i-^q — ± g as follows 



11/ (my ± g,g) 
liRirriT ± g,g) 




. , 0(9 - ?o) - 0(-g - go) T 

4g I Tx^ q 



In 



+ 1 



+ 



TTln 



where we have used eqs.( ^.3|) and (|5.4|) . Thus as anticipated by the discussion of the effective 
static scattering amplitude in the near critical case ( |5.1| ) we see that the real part of the 
polarization is indeed similar to the one in the critical case for the relevant loop momenta 
up to corrections of order mj'/q -C mT/XT in the region of semisoft loop momenta. 

Therefore neglecting terms of order mT/XT which are negligible in the region of interest, 
we find for the spectral densities expressions similar to these in the critical case : 



4g 

p{mT + ^g, 1) = ^ 



L{q) 



XT ^ ^ 



where we have introduced 



+ L^{q) 



L{q) = log 



p{mT 



2q 



11 



Aq 
AT 



Aqn 
AT 



+ mq) 



e rriT 



An analysis similar to that in the critical case reveals that soft loop momenta q <^ XT are 
effectively screened by the renormalization of the scattering amplitude in the near critical 
region, leading to a contribution of order mrp to the damping rate. Semisoft and hard loop 
momenta q > XT are not screened by the resummation of the scattering amplitude and 
determine the leading contributions to the damping rate. 

As argued above, the dominat loop momenta in the integral (|5.5| ) are of order q > q* = 
m 2 3> 7— ^ rriT hence we can approximate 



An 



n 



(ujq) — n{ujq — m-jp) — n{ijjq + m-r) — niuq) ~ n\q) rriT 



iOq ~ 



q 
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allowing an analytical estimate the damping rate from the approximate expression 



To{mT,T) 



NT 



dq 
2^ 



L{q) 



XT ^ ^ 



4g7r 
AT 



+ L{qy 



(5.6) 



The resulting integrals are infrared finite but having recognized that the leading contri- 
bution arises from the semisoft loop momenta q > XT we have introduced an explicit infrared 
cutoff a = C rriT with C ^ 1. Since the integral is dominated by semisoft and hard loop 
momenta g* ^ AT we can approximate further ^ + 1 ~ ^ whence the two contributions to 
the damping rate coincide. Moreover, the dependence in a is negligible in the critical limit. 
This is confirmed by our numerical analysis which uses the exact expression (|5.5| ), the results 
of which are displayed in figures 8 and 9. The integrals in ( |5.(j| ) again reveal a crossover scale 
q* at which 4:7rq*/XT ~ L{q*). For q <^ q* the logarithmic term L'^{q) dominates in the 
denominators and for q ^ q* the term (47rg*/ AT)^ dominates and the integrand falls off just 
as in the perturbative two loops case. 

We now distinguish between the following two possibilities: 

• AT/ (AiT)e~'^'^^^ <^ rriT -C AT/ (47r), to which we refer as the semicritical regime. In this 
case g* -C T is soft and the classical approximation to the Bose-Einstein distribution 
functions applies. Furthermore, we can expand in inverse powers of the logarithm 
log A3 (my) = log and a straightforward analysis along the lines presented for the 
critical case reveals that the damping rate is approximatively constant and given by 



ro(mT,T) 



AT 

27tN 



1 + 



In A3(mj') 



This is the same result as in the case T = T^, k 0, Eq. (|4.15| ). This result is of 



course expected, in the semisoft region of loop momentum T ^ q ^ XT and for 
AT/(47r)e-^'^/^ < mr < AT/(47r) the screening of the scattering amplitude is ineffec- 
tive and the term (g/AT) dominates over the logarithms, therefore the dependence on 
the mass is negligible in this region. 



mr < AT/(47r) e 



-47r/A 



which we refer to as the ultracritical regime. In this case 



g* ^ T is hard and the classical approximation breaks down and the full Bose-Einstein 
occupation factors must be kept. In this regime the logarithm gives the dominant 
contribution in the denominators and the asymptotic damping rate is given by 



ro(mT,T) 



4 



vriVT 



dq (f 



n\q) 



47r 



T 



e rax 



3iV log(T/mT) 



O 



In^ 

rriT 



(5.7) 



which is exactly the same result as in equation ( |4.17| ) with the momentum scale k 
replaced by the thermal mass my. It must be noticed that in this regime there is no 
dependence on the coupling A. 
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Thus we conclude that the relaxation rate near the critical point |T — ^ AT^ and 
k = has the same features as the critical rate for 7^ and T = Tc, provided we exchange 
the infrared scales rriT and k. 

Discussion of the results: 

The two loops calculation in perturbation theory at the critical point revealed the im- 
portance of the different scales of loop momentum. The loop integrals are dominated by the 
contribution of the soft momentum scales q <^ AT. The contribution from semisoft loop mo- 
menta T ^ q ^ XT is subdominant by a factor (k/XT) ln[AT/A;] ^ 1 in the long-wavelength 
limit k <^ AT and the contribution from hard loop momentum modes g > T is suppressed 
even further in the weak coupling limit by an extra power of the coupling A. 

The large N limit leads to a non-perturbative resummation and results in an infrared 
renormalization of the static scattering amplitude as a consequence of dimensional reduction 
and crossover to an effective three dimensional theory for momenta q ^ AT. 

The effective three dimensional coupling that emerges from this analysis of the static 
scattering amplitude at the critical point is X^^effil) = AT/(4g + AT) which is driven to the 
Wilson-Fisher three dimensional fixed point as g ^ 0. Thus soft loop momenta q <C AT 
are effectively screened by this infrared renormalization of the coupling but semisoft loop 
momenta q > XT are coupled with the three dimensional coupling A3(g) = XT/Aq and 
infrared screening is ineffective for these. The importance of this effective coupling for 
the damping rate can be understood intuitively from figures 2, 3 and 4: the resummation of 
bubbles that leads to the effective scattering amplitude also renormalizes the spectral density 
that determines the self-energy, as shown in figure 4. 

This is precisely the most important mechanism that leads to our results in the large N 
limit. Whereas the lowest order perturbative calculation was dominated by the contribution 
of soft loop momenta q -C AT these are effectively screened by the infrared renormalization 
of the coupling which is near the three dimensional fixed point. The dominant contribution 
now arises from the semisoft q ^ AT and hard q > T loop momentum. In the perturbative 
computation at two loops these scales provided subleading contributions of order AT and 
A^T respectively to the damping rate. 

Clearly the resummation via the large N approximation incorporates the screening of the 
soft loop momentum scales but also reveals the emergence of an ultrasoft scale kus ~ ^ . 

At the critical temperature for k 3> k^s or for homogeneous fluctuations near the critical 
point Tc ^ rriT ^ kus the damping rate is determined by the contribution of the classical, 
semisoft loop momentum scales T ^ q ^ XT, the soft scales AT ^ q being screened by the 
infrared renormalization of the coupling and the crossover to a three dimensional effective 
theory. For long-wavelength fluctuations of the order parameter at the critical temperature 
with k < kus or for homogeneous fluctuations near criticality in the ultracritical region 
fTiT < kus the classical approximation breaks down and the damping rate is completely 
determined by the hard loop momenta q > T. Critical slowing down, i.e. the vanishing of 
the relaxation rate as A; ^ or for /c = as T — > T^ only emerges in this ultrasoft limit as 
shown by equations ( [4. 171 ) and (|5.7]). 



Thus, whereas the large N expansion has provided a consistent resummation and the 
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important ingredient of screening of the couplings for the soft loop momentum modes and 
leads to critical slowing down of long wavelength fluctuations important limitations of the 
results obtained here remain. As we argued in the beginning sections a quasiparticle interpre- 
tation of the long-wavelength collective excitations of the order parameter requires that the 
resonance parameter T{k,T)/u!p{k) <^ 1 with ujp{k) being the position of the quasiparticle 
pole or effectively the microscopic time scale of oscillations of these fluctuations. To leading 
order in the large N limit LOp{k) = ^Jk"^ + mj. in the calculation of the damping rate and the 
resonance parameter. Although critical slowing down emerges from the large limit, we 
see that our results to this order indicate that V{k,T) /ujp{k) ^ 1 for k,mT — > 0. 

There are several possible alternatives: either the quasiparticle picture is not appropriate 
to describe the collective fluctuations of the order parameter at or near the critical point 
or further resummations and or other contributions must be taken into account to obtain a 
description of critical slowing down of collective fluctuations that can be understood within 
a quasiparticle picture. In particular an assessment of i) vertex corrections, and ii) wave 
function renormalization must be pursued which, however, are beyond the leading order in 
the large A^ studied here and thus outside the scope and goals of this article. We are currently 
studying these contributions and expect to report our conclusions in a forthcoming article. 



VI. CONCLUSIONS AND FURTHER QUESTIONS 

In this article we have began the program of studying transport and relaxation at and 
near the critical point in second order phase transitions. The focus here is to provide a 
systematic study of critical slowing down from first principles in a phenomenologically mo- 
tivated quantum field theory. The ultimate goal of this program is to assess the potential 
experimental signatures associated with the critical slowing down of long wavelength fluctu- 
ations at the chiral phase transition. Obtaining a robust understanding of such phenomena 
will have important implications in the QGP and or chiral phase transitions in early universe 
cosmology and ultrarelativistic heavy ion collisions. 

Our study reveals novel phenomena that require a non-perturbative framework for their 
consistent and sistematic treatment which is different from the hard thermal loop resumma- 
tion program used in gauge theories. 

Whereas critical slowing down has been studied thoroughly in classical critical phenom- 
ena and these results were used for preliminary estimates of the correlation length at 
freezeout in heavy ion collisions we are aware of only one prior attempt to study 
critical slowing down in a full relativistic quantum field theory and a similar recent analysis 



In references P9| , |30[| the Wilson renormalization group was used to explore the relax- 



ation of the k = mode of the order parameter slightly away from criticality working at 



A^ = 1. The final results of [p9| , p0| are that for t = iT — T^jT^. approaching the critical limit 
the damping rate for homogenous configurations vanishes as logt with v ~ 0.5 — 0.6. The 
analysis of p9| , |30[] , relies on a truncation of the exact renormalization group equations and 



their numerical evolution. In our opinion there are two very important limitations in this 



approach: i) in refs. |2^j3^ the absorptive parts were treated in a rather simplistic manner 
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and associated with the scattering vertex rather than the self-energy, but more importantly, 
this simplified treatment does not include consistently the Landau damping and multiple 
particle thresholds that are the important ingredients in a consistent and sistematic descrip- 
tion of damping and relaxation, ii) the method relies on the numerical evolution of a set 
of equations that had been truncated without a clear control or understanding of the errors 
induced by this truncation. In particular, it is not clear if the result r(t) ~ f^logt is stable 
with respect to different truncations. 

In our study we have systematically focused on the important aspects associated with 
Landau damping, many-particle threshold effects and a consistent study of real-time phenom- 
ena at finite temperature. As is evident in our study of absorptive parts of the self-energy 
in sections II and III a simplified treatment that does not include consistently these can 
hardly reveal the rich hierarchy of scales and the different physics associated with these: 
the soft scale ks — AT/(47r) and the ultrasoft scale kus — ATexp[— 4'7r/(AT)]. Consistently 
to next to leading order in the large N limit we see that slowing down of relaxation of 
long-wavelength fluctuations only begins to emerge at the ultra-soft scale and in contrast to 
the results obtained in [^,^], we obtain T ~ T/[N\ogt]. Thus, although there is agree- 
ment on the statement that the relaxation rate vanishes at criticality the consistent large 
resummation leads to a very different behavior of the relaxation rate. 

Our main results can be summarized as follows: a consistent treatment of critical slow- 
ing down and of transport phenomena at or near a critical point requires a non-perturbative 
framework to resum the contributions from soft loop momenta which is different from the 
hard thermal loop program of abelian and non-abelian gauge plasmas. A perturbative two 
loop calculation reveals clearly the emergence of a hierarchy of loop momentum scales from 
hard g > T to semisoft T ^ g ^ AT and soft AT ^ q, for weak coupling A ^ 1 the 
scales in this hierarchy are widely separated and the semisoft and soft scales are classical. 
Recognizing the shortcomings of a perturbative treatment for long-wavelength fluctuations, 
we implemented a non-perturbative resummation via the large A^ limit to next to leading 
order. The large A^ limit provides a consistent non-perturbative framework for resummation 
of infrared contributions. It clearly displays the infrared renormalization of the scattering 
amplitude in the static limit at or near the critical point and the crossover to three dimen- 
sional physics for soft loop momentum. The resummation of the scattering amplitude leads 
to an effective three dimensional coupling that interpolates between the bare coupling for 
loop momenta q 3> AT and the three dimensional Wilson-Fisher flxed point for q <C AT. 

The infrared renormalization of the effective coupling screens the contribution from soft 
loop momentum to the self-energy and the relaxation rate, which is now dominated by the 
contribution of semisoft and hard loop momenta T ^ q ^ AT. Furthermore a new ultrasoft 
(in the weak coupling limit) non-perturbative scale emerges kus — ^ that signals the 

breakdown of the classical approximation and the dominance of hard loop momentum modes. 

For k , niT ^ kus the damping rate is dominated by the classical semisoft scales and given 
by T{k, T) = whereas for k < kus the hard loop momenta region dominates and leads to 
the damping rate T(k,T) = ^ criticality or TQ{mT,T) = ^^^^^t near criticality for 

k my 
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homogeneous fluctuations, which reveal the slowing down of relaxation of critical ultrasoft 
fluctuations with a damping rate that is independent of the coupling. 

As discussed above, however, these results and those found in p9| , |30| seem to indicate a 
breadown of the quasiparticle picture of collective excitations of the order parameter because 
the resonance parameter r{k,T)/ujp{k) ^ 1 in the long-wavelength limit and the excitation 
decays on time scales much shorter than the natural oscillation time uj~^{k). At this stage 
it is not clear if this feature is a true physical manifestation of relaxation of collective ex- 
citations at or near the critical point or that further resummation and other contributions 
that are beyond the leading order in the large N must be accounted for. We are currently 
studying this possibility by introducing the renormalization group at finite temperature and 
analyzing in detail the contribution from vertex and wave function renormalizations and 
expect to report on further understanding on these issues in a forthcoming article. At this 
stage our study has revealed a wealth of new phenomena and a hierarchy of scales which will 
require a deeper understanding for a complete and consistent treatment of transport and 
eventually hydrodynamics near or at the critical point. Only a thorough understanding of 
these phenomena can lead to an unambiguous assessment of the phenomenological implica- 
tions of critical fluctuations either in the formation of cosmological relics in the early universe 
or in experimental observables in ultrarelativistic heavy ion collisions thus motivating and 
justifying their study. 
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APPENDIX A: EQUATIONS OF MOTION AND SPECTRAL DENSITIES IN 

THE LARGE LIMIT 

We study the relaxation of the order parameter via the real time description of non- 
equilibrium quantum field theory. This formulation requires the time evolved density matrix 
and is cast in terms of a path integral along a contour in complex time, a forward branch 
corresponds to the forward time evolution of the density matrix via the unitary time evolution 
operator and the backward branch represents the inverse unitary time evolution that post- 
multiplies the density matrix. Consequently there are four propagators: corresponding to 
fields on either branch. For a more complete description of this formulation the reader is 
referred to and references therein. The main ingredient in this program are the free field 
Wightmann and Green's functions for the bosonic field In terms of the spatial Fourier 
transform of the bosonic fields <I>^ these are given by 
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(Al) 



^n,S^)^~-k,^^'))o = -^ Kb Gz-{t,t') , 
Gt+(t,t') = G|(t,t') 0(t - t') + Glit,t') e{t' - t) , 

G-f{t,t') = Gl{tX) e{t' -t) + GS{t,t') &{t - t') , 

Gi-{t,t') = Gl{t,t') , 

G.^t,t') = Gl{t,t') , 

GIM = G^^{t\t) . 

(A2) 

where {A{t)B{t')) = Tr[y4(t)i?(t')p(0)] denotes the expectation value of Heisenberg field 
operators with respect to the initial normalized density density matrix which is taken to 
describe a thermal state and the subscript o refers to free fields. It is clear that these real 
time propagators satisfy the identity: 

Gt+{t,t') + G'fit^t') - Gt-it,t') - G-^^{t,t') = . 

The retarded and advanced propagators are defined as 

GA,dt^t') = Gi^tX) - G-^^{t,t') = [Glit,t') - G|(t,t')] e(t' - t) , 
Where for the cases under consideration with fields in thermal equilibrium 

Glit, t') = 4 [[1 + e-^(*-*') + n^, e^'^-'^ , 



G^t^t') 



I 

2k 

i 

2k 



-ik(t-t') 



ik{t-t') 



(A3) 
(A4) 



exp(/5co'^ 



From the Lagrangian density in terms of the auxiliary fields given by ( [4.1| ) it is straightfor- 
ward to find the free field real time correlation functions for the auxiliary fields. In terms of 
the spatial Fourier transform of the auxiliary field Xk these are given by 
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{x-^it)x:^,it'))o = -i5{t-t') 

Figures 10(a) and 10(b) depict the series of Feynman diagrams for the full (x^X^) propagator 
and for the /u// plus-plus component of the propagator of the composite field ($)^ = i.e. 
^($+)2(<|+)2^_ Figures 11(a), 11(b) and 12(a), 12(b) depict similar relations for the {x^X~) 
and (x X ) propagators. Thus using the free field propagators for the auxiliary field given 
above we find the following relations to all orders for the full propagators 



ixtmiiit')) = i5{t - t') + A($t'^(t)$3(t')) (A5) 
{xtmZ^S^')) = A($t'^(t)<f If (t')) (A6) 

(x,^(t)x:g(t')) = -^5{t - 1') + H^-{"m~jlit')) (A7) 



with the definition 

^f^it) = J d^x e'^-^ $±(x,t) ■ $±(x,t) 

The correlation functions of the bilinear composite operator can be written in terms of 
spectral densities in the following manner 

<^>!|(t')) = Jdu; [p>.^.{uj; k) e{t - t') + p<,,M. k) e{t' - t)] ^'^^-^'^ (A8) 
{^t/{t)^-J{t')) = Jduj p< ^.(o.; k) e-^-(*-*') (A9) 

Familiar manipulations introducing a complete set of energy eigenstates in the trace lead to 
the KMS condition 

p<.^.{u;;k)=e^^ p>,^,iuj;k) (AlO) 

The main reason for presenting these formal steps is that the auxiliary field itself does not 
have a KMS relationship for its spectral functions because it is not a canonical field but 
a Lagrange multiplier. However the relations (^-^) which hold to all orders relate the 
correlators of the auxiliary field to those of the bilinear composite operator for which the 
spectral functions associated with their correlators do obey the KMS condition. 

Writing the retarded correlator for the auxiliary fields as a spectral representation 

- = 'J^ e-^^"^'-''\ 



using the spectral representations ( |A8| )-( [A9D and the representation for Q{t — t') given by 
( p.6|) we obtain the relation between the spectral representation for the retarded correlator 
of the auxiliary field and that for the bilinear composite in the following form 
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Px(?o, k) = l + X duj p^2^2iuj; k) 



1 - e"^'^ 
qo- uj + ie 



(Aii: 



where we have used the KMS condition ( |A10 ). The next step of the program is to obtain 
the spectral density Pxilo, k) to leading order in the large N limit. This is achieved through 
linear response analysis for the expectation value of the auxiliary field. 



1. Linear response for the auxiliary field 

The real time expectation value of the auxiliary field is obtained by coupling an external 
source to the auxiliary field in the original Lagrangian ( [4.1|) C ^ C + J^x with the same 
external source for the two time branches. Assuming the addition of counterterms in the 
Lagrangian to ensure that the expectation value of the auxiliary field vanishes for vanishing 
external source we have 



Introducing the Fourier transforms 

m = / ^ Kqo, k) e-^^o* ; J^^^,{t) = J^ J{qo, k) e'^^"* 

we find 

— * — * 

S{k, go) = - Jxik,qo) py^{k,qo) 

We now use the tadpole method [^] to obtain the equation of motion for the expectation 
value of the auxiliary field in leading order in the large limit, thereby obtaining an explicit 
expression for p^i^qo, k) to this order. The implementation of the tadpole method begins by 
shifting the auxiliary field 

= + ; {x{x,t))=Q 

and requiring that {x{^i t)) = to all orders in perturbation theory. A counterterm is added 
to the Lagrangian to cancel the tadpole contributions so as to make the expectation value 
of the auxiliary field to vanish in the absence of the source term thus allowing to extract the 
spectral density straightforwardly. 

To leading order in the large A^ limit we obtain the equation of motion (after the cancel- 
lation of the tadpole term) to be given by 

6{x, t) + j d^x' dt' n^(f -x\t- f) 5{x\ t) = -J-x{x, t) 



with the retarded polarization given by 

. A 



Ur{x-x',t-t')=2t^Y.[{^tiS,mtiS',t)){^^{x,mtiS',t))- 

a,b 
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In terms of the spatial Fourier transform the equation of motion becomes 

/oo 
dt' U,At-t') 5^,{t') = -J^,{t) 
'OO 

and the retarded polarization kernel simplifies to 



= 4A 



(27r)3 
(Pq 



(g+|fc + g1)(t-0 



+ in—n._^^) sin [{\k + q\-q){t- t')] } 0(t - t') 



using the representation of the theta function given by ( |3.(j| ) we find the time-Fourier repre- 
sentation of the retarded polarization to be given by 



-iqoit-t') 



The Fourier transform of the polarization is now written as a dispersion integral in terms of 
the spectral density as 

TT J qo — uj + te 



where 



Uj{q,uj) = 2\tx J 



{2tt)^ Ap\p + q\ 



{[l+n^+j^+Uj;] [5{uj- \p + q\-p)-5{uj+ \p + q\+p)] 

+ [np - rig^+p] [5iuj - \p + g] + p) - +\p + q\- p)]} (A12) 



Finally, in terms of the time Fourier transform the equation of motion for the expectation 
value of the auxiliary field is given by 

^(^0, <?) [1 + n(go, q)] = -JxiQo, Q) 
and we can read off the propagators for the auxiliary field in Fourier space 



(A13) 



The series of diagrams that are being summed leading to the propagator for the auxiliary 
field is shown in fig. (2). 

We now have all of the elements necessary to obtain p^2^2{qo,q), p^2^ 2(go, q), writing 
n(Q'o, q) = ^niqo, q) +i n/(go, q) and comparing the imaginary parts of( |All| ) and (|A13|) and 
using the KMS condition ( |A10| ) we finally find 
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X > / N 1 n7(go,g) [1 + n(go)] ..^.^ 

^ Ml + n«ta.,)]^ + nK*.,) 

X < / A 1 n/(go,g) n(go) 

' " = ^ |n-n«te..)|- + n;fao,.) 

^(^o) 



— 1 



We postpone the evaluation of n(go,9) to Appendix B and now focus on obtaining the 
resummed self energy for the order parameter. 



2. Equation of motion for the order parameter 

We now obtain the equation of motion for the order parameter to 0{1/N) in the linearized 
approximation again via the tadpole method and recognize the self-energy to this order. To 
this effect we write the field as in iU-'^) with 



($^(f,t)) = y.(f,t)5,i ; (r/^(f,t)) = 

where we chose the particular direction "1" by choosing explicitly the external source in 
as .r{x,t) = J{x,t) 5*'^ to give the field an expectation value solely in this direction. 
The equation of motion for tf{x, t) is obtained by imposing that {ri^{x, t)) = consistently in 
the perturbative expansion. In terms of the spatial Fourier transform of the order parameter 
ip{t) we find 

/oo 
J^retAt - i') Vk{t') dt' = Jk{t) 
-oo 

where Jk{t) is the external source that generates the initial value problem and Mf^^{T) ~ 
t)^) ~ 0{N) is the tadpole contribution which is the leading order in the large N 
limit. The 0{1/N) contribution to the self-energy is calculated in terms of the auxiliary 
field and is given by 



with {xtif)x~A^')) s-iid {Xffi't)x^n'it')) the /u// propagators up to 0{1/N) given by (|A5| - [A6D 



in terms of the spectral representations given by with the spectral densities given 



in terms of the self-energy of the auxiliary field by (|A14| - [AT5| ). The contribution to the 
propagator of the scalar field up to order 0{1/N) is depicted in fig. 4. The contribution to 
the auxiliary field propagators from the delta functions ±i6{t — t') gives a local tadpole which 
is cancelled along with the leading order 0{1) tadpole contribution by the counterterm to set 
the theory at the critical point up to this order in the large N expansion. Using the spectral 
representation for the propagators of the auxiliary field and the free field propagators for the 
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bosonic fields given by and after some straiglitforward algebra using the relation 

1 + n{—qo) = —n{qo) we finally obtain 

/oo 
p{lu, k) sin [uj{t — t')] du! 
-oo 

with p{u!, k) given by eqs.( ^^ -( [4.8|) . 



APPENDIX B: THE RETARDED POLARIZATION OF THE AUXILIARY FIELD 



The spectral density (|A12|) is the same as ( p. 51) up to the factor {N + 2)/N, a rela- 
tively straightforward calculation with the Bose-Einstein distribution functions for massless 
particles then leads to 



A I 2T 

n7(go, Q) = ^ { 0(ko| - q) sign(go) + — In 
Svr q 





ko+<7l 1 


— e 


2T 




\<io-i\ 


— e~ 


2T 




1). We 



:bi] 



in the finite temperature contribution to both the real and imaginary part therefore we 
only consider the second term in It proves convenient to write the polarization as a 

dispersion relation 



71 J-oo uj — qo + iO 

and to analytically continue + iqo = s. Using the fact that Ili{uj, q) is an odd function of 
UJ we obtain the dispersion relation 



n(s,g) 



AT /-"^ , UJ 
duj 



2'K^q Jo uj^ + 



In 



1 — e 2T 
1 — e 2T 



A 



In 



+ 



:b21 



where p? is a subtraction point. 

We compute this integral using the sine-Fourier transform as follows. The integrand of 



eq. (p^) is the product of two odd functions of uj : 

/i(^) = 2 /2(^) = In 

UJ"^ + S'' 

We can then apply the Plancherel formula 



1 — e 2T 



1 — e 2T 



c^'^ fii^) f2{^) = dx fi{x) f2{x) 

JO 

where and /2(x) are the sine- Fourier transforms of fi{uj) and f2{uj), respectively. That 
is, 



IT Jo 



duj fi{uj) sinujx 
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where i = 1,2. We find 



and 



We have now that 



f2{x) 



' 2 sin qx 
2T 



2' 



[2'kT X coth(27rTa;) - 1] 



n(s, g) = 7 — — — / dx — — singx Y2'kT x coth(27rTs) — 1] 
{2'kY q Jo x^ 



(2vr)2 q 

It is convenient to split this integral into two terms, 



2 , ^2 



n(s,g) 



A 



00 g~ 
c/x — 



sm qx 



{2'Kfq \Jq 

roc 

+ / dx [2nT coth(27rTx) singa; — 1] 

JO X 



A , q^ + s"^ 
— ^ In — 



We carried out the integration explicitly with the result p7 
A 



n(s,g) 



(2vr) 



-1 + ln 



s g mT 

, - arctg- In 

4tcT q s q 



Y(il±l\ r fl + i 



s+q 



r(fi)r(i + 



is—q 
AmT 



where we used Malmsten formula for the Gamma functions p7 |. 
Back in real frequencies we have, 

n(go ± iO, q) = UR{qo, q) ± i n/(go, g) 

where n7(go, g) is given by eg. (piP and 

A {tx^T 



nR(go, g) 



(2vr)2 

+ ^ln 
2g 



— [0{q-qo)-e{-q-q,)] + \n£^ 



Q-Qo 



2tiT^ , 
H 1mm 



r 1 



AixiT 



r 1 + 



9 + 90 
47riT 



The limit T/g ^> 1 can be taken in a straightforward manner and we obtain the high 
temperature limit of the polarization to be given by 



T-r/ • \ « AT , 
n(go, q + ie) = In 



+ 



Anq 
A 

W) 



qo + ie + q 
qo + ie-q 



AttT 2q \go + «e-g. 



+ 27 



+ o{- 



where 7 is the Euler-Mascheroni constant. 
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FIGURES 




FIG. 1. The function F[x] vs. x = q/k. 
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FIG. 3. Two particle s-channel scattering amplitude in the static and large limit, q, qo = 
are the transferred momentum and frequency carried by the propagator of the auxiliary field. 
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FIG. 4. The self-energy of the scalar field at order 
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FIG. 6. NT{k,T)/T vs. hi[T/k] for A = 1.0. We also plot the classical value A/[27r] and 
the asymptotic ultrasoft behavior NT{k,T)/T = 317^7^- For this intermediate coupling regime 
the classical approximation provides a qualitative a estimate whereas the ultrasoft regime will be 
reached for k/T ^ 3. 10"''. 
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FIG. 7. r(A:,r)/r vs. ln[r/A:] for A = 0.2. We also plot the classical value A/[27r] and the 
asymptotic ultrasoft behavior T{k^T)/T = In the small coupling regime the classical 

approximation describes very well the behaviour of the relaxation rate. The ultrasoft regime will 
only be reached for extremely small momenta k/T <C 8. 10~^^ 
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FIG. 8. Damping rate N ro(mr,r)/T vs. ln[T/mT] for homogenous configurations 
(A; = 0, T ^ Tc) compared with the classical approximation N Tci{mT,T)/T = A/(27r) and 
the asymptotic expression in the medium coupling regime A = 1.0. 
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FIG. 9. Damping rate N ro(mr,r)/T vs. ln[T/mT] for homogenous configurations 
(A; = 0, T ^ Tc) compared with the classical approximation N Tci{mT,T)/T = A/(27r) and 
the asymptotic expression in the small coupling regime A = 0.2. 
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(b) 

FIG. 10. (a): Feynman diagrams for the correlator {x'^ {x,t)x^ (x' ,t')) , (b):Feynman diagrams 
for the correlator ((^+(a;, t))2($+(f', f ))^). 
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(b) 

FIG. 11. (a): Feynman diagrams for the correlator i'))) 0^)'- Feynman diagrams 

for the correlator t'))^). 
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(b) 

FIG. 12. (a): Feynman diagrams for the correlator i'))) 0^)'- Feynman diagrams 

for the correlator f 
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